AV  fcGO  34  & 


CONTRACT  REPORT  NO.  1-170 


A  COMPARISON  OF  PLANE  AND  SPHERICAL 
TRANSIENT  VOIGT  WAVES  WITH  EXPLOSION 
GENERATED  WAVES  IN  ROCK  MASSES 

G.  B.  Clark 
G.  B.  Rupert 
J.  E.  Jamison 


April  1967 
Sponsored  by 

Defense  Atomic  Support  Agency 
NWER  Subtask  I3.I9IC 

Conducted  for 

U.  S.  Army  Engineer  Waterways  Experiment  Station 
CORPS  OF  ENGINEERS 

Vicksburg,  Mississippi 


Under 

Contract  No.  DA-22-079-eng-464 
by 

George  B.  Clark  and  Associates 
Rolla,  Missouri 

VICKSSUSO.  MISS. 


THIS  DOCUMENT  NAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE 
AND  SALE;  ITS  DISTRIBUTION  IS  UNLIMITED 


ABSTRACT 


Plane  and  spherical  waves  in  a  Voigt  medium  were 
investigated  to  compare  the  calculated  wave  forms  with 
observed  waves  generated  by  large  contained  HE  and  NE 
explosions.  Interest  is  centered  on  wave  forms  in  what 
is  usually  considered  to  be  the  elastic  region  around 
an  explosion.  Plane  waves  do  not  apply  at  this  distance 
because  of  geometry. 

The  plane  Voigt  wave  equation  has  been  previously 
solved  for  particle  velocity,  stress  and  strain  for  a 
unit  impulse  forcing  function.  However,  solutions  for 
the  displacement  for  a  unit  impulse  and  for  the  four 
wave  parameters  for  a  unit  step  and  a  decay  exponential 
involve  multipliers  in  the  operational  form  for  which  no 
transform  pairs  have  been  published.  A  method  of  solu¬ 
tion  is  presented  which  utilizes  a  Heaviside  expansion 
of  the  multipliers  in  the  transform  plane  which  results 
in  products  of  two  infinite  series  which  may  be  inverted 
term  by  term.  These  may  be  further  resolved  as  single 
series  with  polynomial  coefficients  for  purposes  of  com¬ 
putation. 

A  similar  method  of  solution  of  Voigt  spherical 
waves  was  found  for  unit  impulse,  unit  step  and  decay 
exponential  forcing  functions  for  displacement,  parti'le 
velocity,  strain  and  radial  stress.  Appropriate  recur 
sion  formulas  make  them  readily  adaptable  to  computer 
evaluation.  Oscillations  occur  for  a  spherical  wave 
whereas  for  a  plane  wave  they  do  not. 

Calculations  were  performed  for  particle  velocity 
for  three  values  of  the  Voigt  viscoelastic  parameter  w0 
and  comparisons  made  with  pulse  forms  for  waves  in  granite, 
tuff  and  salt. 
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A  Comparison  of 

Plane  and  Spherical  Transient  Voigt  Waves 
With  Explosion  Generated  Waves  in  Rock  Masses 


Introduction 


The  Voigt  model  yields  responses  to  transient  forc¬ 
ing  functions  which  are  different  from  the  standard  elas¬ 
tic  model,  and  earlier  (Ref.  1)  mathematical  analyses  of 
plane  waves  in  Voigt  materials  indicated  pulse  lengthen¬ 
ing  and  attenuation  somewhat  similar  to  that  observed  in 
waves  generated  by  large  explosions  in  rock. 

While  the  theory  of  elasticity  provides  a  good  first 
approximation  to  many  problems  in  rock  mechanics,  in  the 
solution  to  close-in  wave  mechanics  in  rock  it  does  not 
give  an  accurate  descripticn  of  the  behavior  of  peak 
stress,  strain,  displacement,  velocity  or  pulse  length. 
The  Voigt  viscoelastic  model  is  among  the  many  physical 
analogs  which  have  been  investigated  in  an  attempt  to 
find  valid  mathematical  representations  of  wave  para¬ 
meters  in  rock.  The  only  solutions  for  a  plane  Voigt 
wave  in  the  literature  appears  to  be  those  of  Collins 
(Ref.  1)  and  Hanin  (Ref.  2)  in  which  they  solved  for  the 
particle  velocity  in  response  to  a  unit  impulse. 

Lee  (Ref.  3)  solved  the  spherical  Voigt  equation 
for  stress  for  a  continuous  harmonic  forcing  function. 
However,  in  solving  for  response  to  transient  loading 
he  introduced  a  "constant  loss  factor”  into  the  trans¬ 
form  which  changed  the  problem  from  a  Voigt  model  to 
essentially  a  solid  friction  model. 

The  objectives  cf  this  investigation  are  to  arrive 
at  and  present  solutions  to  both  the  plane  and  spherical 
wave  equations  for  parameter  of  interest,  i.e.,  particle 
velocity,  strain,  displacement  and  stress,  and  to  compare 
the  calculated  results  with  these  obtained  in  observa¬ 
tions  of  waves  in  rocks  from  large  HI!  and  nuclear  explo¬ 
sions. 

There  have  been  many  approaches  made  to  the  solu¬ 
tion  of  transient  Voigt  waves,  but  apparently  the  first 
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successful  solution  was  made  by  Collins  (Ref.  1)  for 
particle  velocity  for  a  unit  impulse.  Clark  and  Rupert 
(Ref.  4  and  5)  obtained  stress  and  strain  by  Collins' 
method  and  the  displacement  by  numerical  integration. 

The  method  employed  by  Collins  was  to  normalize  and 
shift  the  Laplace  transform  solution  for  the  Voigt 
wave  equation,  expand  the  positive  part  of  a  resulting 
exponential  in  a  series,  and  obtain  a  term  by  term  in¬ 
version  expressed  as  negative  order  parabolic  cylinder 
functions  multiplied  by  a  factorial,  an  exponential  and 
a  power  term. 

Operational  solutions  for  the  displacement  for  a 
unit  impulse  and  for  forcing  functions  such  as  a  unit 
step  or  a  decay  exponential  require  that  the  transform 
solution  also  contain  an  appropriate  multiplier.  N’o 
transform  pairs  for  the  resulting  expressions  exist  in 
published  tables. 

The  transformed  multipliers  may,  however,  be  ex¬ 
pressed  as  Heaviside  expansions,  yielding  3  transformed 
solution  in  the  s-plane  composed  of  the  product  of  two 
infinite  series.  These  may  be  inverted  term  by  term 
into  th"  time  plane,  expressed  as  a  product  of  two 
serie-  in  the  time  variable  and  by  proper  arrangement 
utilized  for  computer  evaluation.  Some  of  the  double 
series  may  be  expressed  as  Cauchy  products,  which  clari¬ 
fies  the  mathematical  operations  and  assists  materially 
in  computer  programming. 

V'oij’t  Model:  The  basic  Voigt  model  for  an  aeolotropic 
soT  iT  was  conceived  for  a  substance  in  which  it  wu. 
assumed  that  the  stress  components  could  oe  expressed 
as  the  sum  of  two  sets  of  terms,  the  first  being  pro¬ 
portional  to  the  strains  and  the  second  proportional  to 
the  rate  of  change  of  the  strains.  {Rif.  t>)  For  an 
aeolotropic  solid  the  stress  equations  arc  of  the  form 
(Ref.  ’) 


1 c  1  r 


dt  1 


*t 


*  Cl  1*  d  1  \  ) 


.  «•  1< 


[ 1) 
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where  a  o„,  -r  and  t  .  are  of  similar  form  with 

y  2  yz ,  xz  x> 

appropriate  c^  and  d^  multiplying  constants.  If  energy 

is  to  be  an  univalued  function  then  c  *=  c  and  d  „  = 

rs  sr  rs 

dgr.  This  reduces  the  number  of  independent  constants 
from  72  to  42.  In  an  isotropic  medium  there  are  only 
three  elastic  constants  and  three  viscoelastic  constants. 
Thus , 

C 12  =  c  1 3  -  *-2 1  =  C31  =  C32  =  X 
dl2  =  d 1 3  =  d2!  =  d  3 1  =  d32  =  X' 


C44  =  C55  =  c66  =  v 

d44  =  =  dg6  =  u' 


(2) 


C11  =  c22  =  c3?  3  X+2p 


dii  =  d22  -  d 3 3  =  X'+2y' 


Thus, the  stresses  for  a  Voigt  solid  become 


where 


°x =  i?)  a  +  2(,j+ij'  y  £x 

°>-  =  (*♦*’  It)  a  +  2(*j+‘j'  It)  cy 

°z  =  (A+A’  It)  a  *  2^*  It)  s 

Tyz  =  (u+u'!t)  V- 

Txz  =  ^+u'lr)  yxz 

Txy  =  (“♦’‘'IT)  V 


A  =  e  *  e  +  e  (dilatation) 

x  y  z 


Thus,  the  wave  equation  may  be  obtained  for  a  Voigt 
solid  by  substituting  (a  +  a'  |-t-)  and  (y  +  y'  |—)  for  v 


1 


I 


Plane  Waves 

The  basic  wave  equation  for  a  plane  Voigt  wave  may 
be  derived  in  different  ways  and  in  terms  of  alternative 
parameters : 

(i*i»)  (4) 

(Do  at''  3x2  c2  3t2 

where  z  represents  either  displacement,  velocity,  stress 
or  strain. 

Where , 

us o  =  viscoelastic  or  transition  frequency  factor 
c  =  bar  velocity  =  E/p 
x  =  distance 
t  =  time 

Collins,  (Ref.  1),  however,  derived  the  operational  form 
of  the  equation  directly  for  particle  velocity  utilizing 
the  force  acceleration  relation 


3o  3v 
3x  =  p  3t 

in  which  particle  velocity  is 


(5) 


(6) 


and  the  relationship  of  strain  to  displacement 


(7) 


together  with  the  following.  The  Voigt  relationship 
between  stress  and  strain  is 


o  =  E(e 


+ 


1_  3e_-> 
wo  St-1 


(8) 
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which  in  operational  form  is 


a  =  E(1+s/6o<,)£ 


(9) 


Equation  (5)  in  operational  form  becomes 


(10) 


where  the  bar  indicates  that  the  Laplace  transform 
has  been  taken  and  (6)  is 

u  =  sv  (11) 


For  a  delta  function  input  the  boundary  condition  at 
t>o  is 

a  =  -  P ' 5 (t)  (12) 

which  is  unity  in  operational  form,  i.e.. 


o  =  -P' 

From  (10)  and  (11) 

3o 

37=  ps  u 

or  differentiating 


(13) 


(14) 


3  2a 
3x2 


ps' 


3u 

3x 


ps' 


(15) 
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0 


which,  with  (9)  becomes 


3a0  _  £_  s2  - 
ax2”  ‘  E  (s/aio+l) 


(16) 


The  solution  to  this,  after  satisfying  the  boundary 
conditions  is 


a  =  exp[-xs/c(s/w0+l) 


(17) 


Differentiating  once  with  respect  to  x  and  using  (10) 
gives  the  velocity  equation 

-  P'  exp[-xs/c(s/u„+l)*]  (18) 

pC  /s/wo+1 


The  notation  of  Collins  and  van  der  Pol  is  used 
herein  for  the  symbol  connecting  a  function  to  its 
transform: 


v(x,t)  v(x,s) 


.19) 


The  transform  expression  in  ''  0  is  then  "normalized" 
by  substituting  ni,s  for  w0,  or  multiplying  s  by  w0  •  This 
is  equivalent  to  dividing  the  time  and  the  velocity  func¬ 
tion  by  wo,  thus 


i_  V(x,t/<D0J  -  exjPL-W^Ls+lj*]  (20) 

“#  'pc  /s+1 


The  shift  theorem  is  then  used  to  place  the  transform 
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solution  in  a  form  which  can  be  inverted  by  use  of  tabu¬ 
lated  transform  pairs,  i.e., 


v(x,t/iD0)e 


P'uo  exp[-X(s-l)/sJ] 

pc  7s 


(21) 


where  X  = 

c 


The  second  half  of  the  exponential  is  then  expanded: 


n=o 


£  1 

n!  772 


(22) 


which  gives 


> 

'■..i 


v(x,t/'wo)et 


-  P'M° 
•  pc 


l 


n=o 


exp(-xTs) 
n!  sn/2+l/2 


(23) 


By  use  of  formula  (9)  page  246  of  Ref . (8) equation  (23) 
may  be  inverted  term  by  term  to  give 


/  *  /  s  t  P  1  (Dp 

v(x,t/u)o)e  *  — 

P  ^ 

°-n  ^ 


n=o 


xn2n/2tn/2-l/2 
n !  A 


exp(-X  2/8t) 


(24) 


or,  rearranging  and  substituting  for  t  and  letting 
(00t  =  T 


V(x,t) 


X"  ?n/2  (n-l)/2  ,  X  x 

n! 


(25) 
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Similarly,  from  equation  (24) 


o(x,t/oj0)et  *=  -  P'w0  exp[-X(s-l)/s4) 


(26) 


or 


o(x,t/w0)e 


t 


P'Uo  l 
n=o 


(27) 


whose  solution  is 


o(x,t) 


-P>»  V  xV(n-l)/2  n/2-1  r  X  > 

/it  ^  n!  -n+1  '■721-' 

n=o 


(28) 


The  displacement  is  the  integral  of  the  velocity  with 
respect  to  time,  which  is  obtained  by  multiplying  equa¬ 
tion  (11)  by  1/s.  This  operation  is  valid  only  for  zero 
displacement  at  t  =  o.  Hence 


u(x,t) 


P'  exp[-xs/c(s/w0+l)  j 
pc  s  /s/a'  +1 


(29) 


Normalizing  and  shifting  gives 


—  u(x,t/(Do)et 
Wo 


P_|_  exp[-X(s-l)/s^ 

pc  Wo  (S-  1)  /s" 


1 30) 


In  this  case  l/(s-l)  must  be  expanded  and  the  final 
operational  solution  is  a  double  series,  i.e., 


1 _  _  l  r  l_ 

s-1  s  “  m 
m=o  s 


which  yields 


t.  r'  “  v  Xn  exp(-x/s) 

u(xtt/w0)e  l  1  fH"  n)2+m+V? 

m=o  m=o  s 


which  may  be  inverted  term  by  term  to  give 

p'  0  r  r  Xn  1n/2+m+l_n/2+m=l/2 

u(x,t)  =  i—r  exp(-T-X2/8T)  \  l  2  T 

m=o  n=o 


l)  ,  ,,  X  , 

-n-2m-2(-^p) 


The  double  series  is  readily  amenable  to  evaluation 
by  computer  methods  for  small  values  of  X  and  T. 

The  expression  for  strain  is  found  by  differentiating 
equation  (29)  with  respect  to  x,  which  gives 


^  P'  exp[-xs/c(s/u)o*l) 

t(s>t)  =  ^7"  isTOii - 
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which  yields 


1  ,  t  ,  t 

—  e(x,— )  e 

iii  o  U)© 


(35) 


Equation  (35)  may  be  inverted  term  by  term  to  give 

e(x,t)  *  — j-p  exp(-T-X2/8t)  I  jjT  2(n+1^V/2D_n  (36) 

pcVn  n=o  S2j 


The  solutions  for  an  exponential  decay  input  can  be 
obtained  by  multiplying  the  transform  for  each  parti¬ 
cular  parameter  before  normalizing  and  shifting  by 
Ej/(s  +  e),  wl  ich  is  the  transform  of  P(t)  =  P0exp(-6t). 
This  is  illustrated  by  the  derivation  of  the  velocity 
equation: 


v(x,t) 


1 

(s*5'y 


,  X  3 

exp  ( ™ 


/s/Wo* 1 


(37) 


Upon  shifting  and  normalizing  this  becomes 


1_ 

(JO 


v(x,— ) 

Wo 


e 


pc 


1 


a,o  (S-(1-^-)  ] 

UTO 


exp[-x(s- 1)/ /s' ] 

f— 

/s 


(38) 


Let 


[s-(l-  ) ]  «  (s-a) 


(39) 
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and  expanding  the  positive  exponential  as  previously 
done,  and  letting 


1 _ 

s-a 


00 


m=o 


(f>- 


The  transform  then  becomes 


1_ 

Ulo 


)et 


Po _ r  ,a.m  r  X^  expj-Xi'T) 

pcujo  s  ^  ^  n!  n/2<  \/2 

m=o  n=o  s 


which  may  be  rewritten 


Po 

P  CU)o 


00  CO 


l  l 


r 

n! 


m  exp  (-X/T) 
a  sn/2+ra+ 1/2 


The  inverse  is  found  term  by  term  to  be 


v (  St) 


P. 


pc 


cxp(-T-X*/8T)  [  [ 


m=o  n=o 


m,n/2»m» l_n/2  +  ia*  1/2 
a  -  I  i’ 


which  can  be  evaluated  by  computer  methods. 
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the  operational  solution  of  equation  (44)  becomes 


r$(r,t)\=  - 


» 

>  'To  expt-sC'r-rol/cCs/ua*!)7] _ 

rus  (s/-),  +  l)  3s2  + _ 4s _ 

c2(s/u0+l)  r<,c(s/w0  +  l)  i 


(47) 


However,  interest  is  centered  '•n  parameters  such  as 
displacement,  particle  velocity,  stress  and  strain. 
The  displacement  is 


=  IF  (48) 


or 


u(r  o  V.  P 1  r.  expj-sX^r.Vc (s^/w.*! ) 1  +  _ s _ 

u(s/ojo*1)  [Bj  1  r2  rc(s/uj.*!): 


(49) 


where 


(B| 


r  3sJ 
cJ(s/w.»l) 


4s 

r.c  (s/u). *  1 )  j 


(50) 


Let  s  -  *8s  (normalize)  and  then  s  -  s-l(shift), 
and  equation  (49)  becomes 


u(r,t/v. )el 


P’r.  cxp( -R>7«R//s)  jJ1_  t  ( s - 1 )  | 

c».’!Aj  r:  rc»s 


vSM 


2) 


where 


[A]  -  [s2  ♦ 


4c 


3row0 


3/2 

s  + 


4c^ 


3r02w02 


-2)s  -  - 


4c 


*37-0  Ulo 


i] 


(52) 


and  R  =  oi0(r-r6)/c  (53) 

The  fraction  involving  the  quadratic  in  s  may  be 
expanded  as  follows: 


1/[A]  =  ~ 
sz 


I 

m=o 


(54) 


where  the  coefficients  are  functions  of  the  coeffic¬ 
ients  of  the  quadratic.  ‘  The  positive  exponential  may 
likewise  be  expanded  to  yield: 


exp  R//s  =  l 
n=o 


1 

n!  (s)n/2 


(55) 


Substituting  in  (51)  and  separating  into  appropriate 
terms  gives  the  operational  solution  in  terms  of  double 
infinite  series: 


u(r,t/ujo)e 


t 


P'r0c2 

ywo 


00  <30 


l  l 


Rn  .  exp(-R/s) 
n!  m  sn/2^m/2+2 


(56) 


[ 


00  00 


l  l 


b 


m 


exp(-R/s)  V  V  Rn  u  exp(-R/s)  , ] 

_n/2+m/2+3/2  1  1  n!  m  n/2+m/2+5/2  J’ 

s'  m=o  n=o  s 


22 


Equation  (56)  may  then  be  inverted  term  by  term  (Ref.  7) 
to  yield  the  solution  in  the  time  plane. 

The  velocity  for  a  unit  impulse  may  be  found  by 
differentiating  the  displacement  with  respect  to  t,  or 
multiplying  (49)  by  s.  The  strain  is  determined  by 
differentiating  the  displacement  equation  (49)  with  re¬ 
spect  to  r,  and  the  stress  by  substituting  into  equation 
(46).  All  of  these  lead  to  solutions  consisting  of  double 
series  similar  to  equation  (56),  which  may  be  inverted 
and  evaluated  by  computer  methods. 

For  example: 

u(r,t)  =  exp(-T-R2/8T) 


Tn/2.m/2+2D  (-=)  ♦  4 

•A  _m_Q  v  r. _ J  I 


.  ec-£=3  ♦  — l  y  y  5lb  sn/2+tn/2+2 

-n-m-5  /2T  rc|  ^  ^  nl  in 

[0=0  n*o 

Tn/2+m/2*3/2D  ,  Rj  ^  —  b  2n/2+m/2+3 

'n-m-4V/2T  m=o  n=o  m 

Tn/2>m/2+3/2D  (  R  ]  j 

-n-m-6  V2T  ) 


Tn/2>m/2+3/2D  (R 

-n-m-6  V2T 


where 

T  =  Wot 

D-n-m-z  =  Weber  functions  of  negative  integral 
order 


The  expansion  of  fractions  has  been  shown  to  be 
valid  for  both  rational  (Ref.  8)  and  fractional  expo¬ 
nents  (Ref.  9)  of  the  denominator  terms  of  the  frac¬ 
tion  (see  below).  Some  of  the  expansions  in  series  con¬ 
verge  slowly  for  large  values  of  dimensionless  distance 
and  dimensionless  time.  The  limits  of  convergence  for 
real  time  are  therefore  dependent  upon  the  value  of  w0, 
as  well  as  the  programming  techniques  employed  and  the 
capacity  of  the  computer.  Several  other  methods  of  ob¬ 
taining  a  function  which  could  be  inverted  were  investi¬ 
gated  but  none  offered  a  means  of  solution. 

In  every  term  in  the  inverted  equations,  there  are 
essentially  five  components  an  exponential,  a  power- 
factorial  term,  2  to  an  exponential  value,  T  to  an  ex¬ 
ponential  value,  and  a  cylinder  function  of  negative  in¬ 
tegral  order.  The  exponential  decreases  rapidly  with  in¬ 
creasing  values  of  T,  and  the  cylinder  functions  also  de¬ 
crease  with  an  increase  in  order  or  an  increase  in  value 
of  the  argument.  The  power- factorial  term  increases 
rapidly  until  n  =  R  and  then  it  decreases.  Cylinder  func¬ 
tions  are  calculated  by  means  of  an  appropriate  recursion 
formula  (Ref.  4). 

The  factor  (2T)n  increases  without  bound  as  n  in¬ 
creases.  The  behavior  of  all  of  these  functions  must  be 
considered  in  programming  inasmuch  as  some  of  the  numbers 
may  become  very  large,  in  excess  of  1010®,  and  some  very 
small,  less  than  lO'lOO.  in  some  cases  the  order  of  mul¬ 
tiplication  becomes  important  so  that  two  small  or  two 
large  numbers  are  not  multiplied  in  succession. 

For  computation  of  the  double  series,  it  is  more 
convenient  to  find  the  multipliers  of  each  successive 
cylinder  function.  These  are  found  to  be  increasing  trun¬ 
cated  series  of  power- factorial  terms.  That  is,  a  double 
jressed  as  a  Cauchy  product: 

2(n+m+z)/2T(n+m+z-l)/2[) 

-n-m-z 

(58) 

}  2(n+z)/2T(n+z-l)/2D 
'  -n-z 
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y  y  b  K 
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«  n  n-k 

I  (  l  k 
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The  infinite  series  all  converge  rapidly  for  relatively 
small  values  of  R  and  T,  but  not  so  rapidly  for  larger 
values.  Specific  values  must  also  be  chosen  for  r,  r0, 
too,  and  c,  which  was  not  necessary  in  the  case  of  plane 
waves . 


A  rational  function  of  the  following  type  may  be 
expanded  (Ref.  9)  in  an  infinite  series,  and  has  a  zero 
or  order  n  at  ■*>,  i.e.. 


y(s)  = 


"n  n"n 
a0s  +ais  +--a 


may  be  expanded  as 


1  al  1  ai-a0a2  . 
y  (s)  =  - 5-  — r-  + - - - — s- 1— 


a°Sn  a°  s 


2  n+1 


So  ' 


n+2 


and  inverting  term  by  term  yields 


tn-i  tn  tn+l 

y(t)  =  Co  orr rr+  ht+  c2  -&njr+- 


(59) 


(60) 


(61) 


However,  the  individual  coefficients  of  the  terms  in 
the  series  rapidly  become  very  cumbersome  to  develop  and 
to  use  in  computation.  For  digital  computer  calculations, 
on  the  other  hand,  a  simple  sequential  loop  procedure  can 
be  employed  to  determine  the  successive  coefficients. 


Let  equation  (60)  be  written  in  the  following  in 
finite  series  form: 


y(s) 


Co  ci  c2  n 

—  +  - -  +  - —  +  --  — * —  + 

n  n+1  n+2  2n 

s  s  s  s 


The  same  process  may  be  carried  out  for  irrational 
fractions  (Ref.  10)  : 


(62) 


y(q)  = 

where 


1 


n  n-1 
a0q  +aiq  +- 


-a 


(63) 
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This  may  likewise  be  expanded  in  the  form 


and  may  also  be  inverted  term  by  term  to  give  a  solu¬ 
tion  in  the  time  plane. 

Thus,  the  inversion  of  the  double  series  resolves 
itself  into  the  inversion  of  a  single  series  with  a 
polynomial  coefficient,  each  term  of  the  series  having 
a  valid  inversion.  The  inverted  terms  are  similarly 
expressed  and  are  in  a  convenient  form  for  computation. 

The  operational  solutions  for  unit  step  and  decay 
exponentials  are  obtained  in  a  similar  manner.  In  these 
cases  the  operational  expression  is  multiplied  into  the 
quadratic  in  s  before  it  is  expanded. 

It  is  notable  that  there  is  an  oscillation  of  all 
parameters  for  the  spherical  wave,  while  the  plane  Voigt 
wave  does  not  exhibit  oscillations  (Ref.  3) .  The  be¬ 
havior  is  somewhat  similar  to  an  elastic  wave,  but  with 
greater  damping  and  a  difference  in  wave  shape.  For  a 
unit  impulse  and  a  decay  exponential  forcing  function 
the  oscillations  are  about  the  zero  axis.  For  a  unit 
step,  however,  the  oscillations  are  about  a  curve  which 
is  parallel  to  the  zero  axis.  In  an  elastic  material, 
a  wa^e  caused  by  a  unit  step  function  also  oscillates 
about  a  line  parallel  to  the  zero  axis. 

For  larger  travel  distances  the  wave  spreads  out 
and  becomes  somewhat  more  symmetrical.  However,  very 
small  disturbances  which  are  characteristic  of  a  "dif¬ 
fusion  model"  are  indicated  before  the  arrival  of  the 
main  wave. 


Computation  Procedures 

A  generalized  flow  diagram  for  the  computer  program  is 
given  in  Fig.  I  .  The  cylinder  functions  are  evaluated  by 
use  of  the  following  recursion: 


(n+l)D_n2(z)  =  -  zD_nl(t)  ♦  D_n(z) 

(65) 

where 

D.M  .  e-'2/4 

(vIM 

and 

d  - 1  ( z)  =7k2/4  erfc  7= 

(67) 

The  complementary  error  function  is  computed  from 

erfc  =  1.0-erf  (68) 

and  the  error  function  may  be  determined  from  several  ap¬ 
propriate  formulas  (Ref.  11),  including  an  infinite  series 
expansion.  For  large  arguments  the  error  function  rapidly 
approaches  a  value  of  1.0  and  the  value  of  erfc  becomes 
very  small,  which  restricts  the  use  of  equation  (68  )  and 
the  upper  limit  of  argument  which  may  be  employed  in  com¬ 
putation  procedures,  i.e.,  above  a  value  of  about  3.0. 

The  argument  of  the  cylinder  function  is 

R//2T  =  w0 (r- r0) /c/2w0t ,  which  means  that  any  combination 

of  the  quantities  in  the  expression  such  as  large  r,  large 
wo  or  very  small  t,  which  yields  values  of  the  argument 
larger  than  3.0  leads  to  computing  difficulties.  Further, 
large  values  of  w0  in  the  expression 

-T-R2/8T  ,  -co0{t-(il^)2/8t}  ...  ,  - 

e  *  1  v  c  '  1  yield  very  small  values  for 

this  function.  The  cylinder  functions  of  high  order  also 
are  verv  small.  On  the  other  hand  the  factor 

(2T)n«(2w0t)n 

increases  very  rapidly  for  large  values  of  w0  to  3000  for 
the  ranges  of  r  and  t  shown  in  Fig. 33  .  Also,  the  terms 
in  the  polynomial  in  R  increase  very  rapidly  for  large  values 
of  wo  (Eq.  58). 

Thus,  the  series  expressions  are  well  adapted  for  small 
values  of  the  parameters  involved  with  the  exception  of  very 
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small  t.  Numerical  integration  of  the  wave  equations  was 
also  considered,  but  this  method  also  becomes  unwieldy  for 
large  values  of  oj0. 

Comparison  of  Voigt  Waves  and  Observed  Waves 

The  important  parameters  of  waves  in  natural  earth 
materials  are  arrival  time,,  rise  time,  fall  time,  rate  of 
attenuation  with  distance  and  oscillating  characteristics. 
In  this  study  the  observed  waves  of  interest  are  radial 
pulses  generated  by  underground  nuclear  explosions,  i.e., 
in  the  HARDHaT,  AIRVENT,  LONG  SHOT,  and  GNOME  events. 

There  are  several  types  of  energy  losses,  dispersion, 
etc.  in  the  transmission  of  dilatation  waves  through  natu¬ 
ral  rock  masses  which  cause  the  waves  to  differ  from  those 
predicted  by  elastic  theory.  Two  immediate  evidences  of 
this  are  greater  attenuation,  and  lengthening  of  the  ini¬ 
tial  pulse  with  distance.  In  many  rocks  the  pulse  length 
has  been  found  to  increase  approximately  with  the  first 
power  of  travel  distance.  As  indicated  earlier  one  of  the 
primary  purposes  of  this  investigation  was  to  ascertain 
whether  the  spherical  Voigt  wave  equation  would  demonstrate 
any  of  the  characteristics  found  in  observed  pulses  created 
by  under  ground  nuclear  explosions. 

Since  there  have  been  only  approximated  values  for 
for  rocks  published  in  the  literature,  a  range  of  values 
was  used  which  would  demonstrate  the  possible  applicability 
of  the  Voigt  equation.  The  values  employed  for  wD  were 
600,  2000,  and  3000,  having  the  dimensions  of  t'*.  Larger 
values  could  not  be  handled  by  the  computer  program,  and 
smaller  values  represent  large  effective  coefficients  of 
viscosity  which  do  not  appear  to  be  applicable.  That  is, 
the  higher  the  viscosity  of  the  dashpot  element,  the  more 
nearly  the  Voigt  model  acts  as  a  single  rigid  member.  This 
results  in  very  high  velocities  of  the  wave  front  which 
velocities  are  not  characteristic  of  rock  or  other  natural 
materials . 

For  lower  viscosities  the  Voigt  wave  forms  do  demon¬ 
strate  some  of  the  important  properties  of  natural  waves, 
i.e.,  attenuation  and  pulse  lengthening.  For  high  viscosi¬ 
ties,  (small  iD0),  they  exhibit  velocities  of  the  wave  front 
which  are,  as  indicated  above,  velocities  which  exceed  the 
velocity  of  both  elastic  and  observed  waves  in  rock. 
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Plane  Waves.  Plane  Voigt  waves  were  investigated  in  the 
process  of  developing  procedures  for  solving  the  transformed 
spherical  wave  equation.  Close  to  the  source  of  disturbance 
these  waves  demonstrate  a  significant  attenuation  of  peak 
velocity  and  strain,  etc.,  whereas  a  plane  elastic  wave  is 
propagated  without  loss  of  amplitude  or  change  in  shape. 

While  the  behavior  of  plane  Voigt  waves  (Figs. 2  to  13)  is 
not  of  direct  importance  in  this  study,  the  principles 
evolved  in  their  solution  were  valuable  aids  in  solving  and 
programming  the  spherical  wave  equations. 

Forcing  Functions.  As  previously  discussed,  three  forcing 
functions  were  employed  in  calculating  wave  forms,  a  unit 
impulse,  a  unit  step  and  a  decay  exponential.  Each  of 
these  functions  makes  the  transformed  solutions  successively 
more  complicated,  and  the  first  two  were  used  primarily  as 
a  means  of  developing  procedures  for  obtaining  inverse  trans¬ 
forms,  working  from  the  less  to  the  more  complex  forms. 

For  simulation  of  a  real  pressure  function,  the  decay  ex¬ 
ponent  is  the  more  realistic  of  the  three,  but  has  too 
sharp  a  rise  time  for  most  field  conditions.  Finally,  solu¬ 
tions  obtained  for  the  unit  impulse,  6(t),  forcing  function 
serves  as  a  basis  for  solutions  for  other  forcing  functions 
by  means  of  convolution. 

Rock  Properties.  Available  data  on  properties  of  rocks  which 
constituted  the  environment  of  the  nuclear  events  discussed 
are  listed  in  Table  I.  *  data  were  available  for  AIRVENT. 
The  data  for  LONG  SHOT  are  limited  to  strength  and  ultra¬ 
sonic  pulse  velocity.  Data  on  salt  and  granite  are  more 
complete.  However,  no  values  for  viscoelastic  properties 
are  available  inasmuch  as  meth  .  have  not  been  devised  for 
measuring  such  properties  of  rock  masses. 

Physical  property  data  and  viscoelastic  wave  analyses 
are  not  complete  enough  to  derive  any  quantitative  rela¬ 
tionships  between  them.  However,  for  rocks  with  lower  modu¬ 
li  and  a  significant  porosity,  the  pulse  forms  are  of  shapes 
similar  to  those  obtained  from  a  Voigt  model.  Also,  the 
Voigt  wave  shows  some  lengthening  with  travel  distance,  but 
is  not  as  great  as  that  found  in  natural  rock. 
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Wave  Velocity,  A  detailed  frequency-wave  velocity  analy¬ 
sis  has  not  been  made,  but  an  inspection  of  the  particle 
velocity  pulses  indicates  that  the  macro-wave  front  veloci¬ 
ty  for  we  *  600  (Fig.  27,  for  example),  2000,  and  3000  is 
greater  than  the  elastic  velocity,  decreasing  with  an  in¬ 
crease  of  u„.  The  actual  first  disturbance  at  the  distances 
investigated  are  extremely  small,  and  they  occur  at  very 
small  times.  That  is,  there  is  no  well  defined  arrival 
time  for  either  plane  or  spherical  Voigt  waves.  This  is 
not  the  case  for  an  elastic  wave,  however,  where  the  arri¬ 
val  time  is  defined  by  the  elastic  velocity,  i.e.,  for 
t'  *  t-(r-r0)/c,  and  no  wave  is  defined  for  i  1  0. 

For  larger  values  of  u0  the  velocity  at  the  wave  front 
decreases,  but  in  each  case  approximately  one- fourth  of  the 
"measurable"  use  time  has  passed  before  the  corresponding 
elastic  pulse  would  have  arrived.  For  very  large  viscosity 
the  Voigt  model  would  indicate  a  small  disturbance  at  •  for 
a  very  small  time.  For  low  viscosities  th ?  behavior  approaches 
that  of  an  elastic  wave. 

The  forcing  function  of  greatest  interest  in  this  study 
is  the  decay  exponential,  Poe'S*.  The  early  arrival  times 
of  waves  for  this  function  are  in  part  due  to  the  abrupt 
(zero)  rise  time  in  pressure.  The  dashpot  in  the  viscoelastic 
model  tends  to  respond  as  a  rigid  segment  to  instantaneous 
change  in  pressure  or  force.  Thus,  in  a  Maxwell  element  the 
first  response  to  a  sharp  pressure  front  is  almost  totally 
in  the  spring.  In  the  Voigt  model  where  the  spring  and  dash- 
pot  are  in  parallel,  however,  one  element  cannot  respond  with¬ 
out  the  other.  Consequently,  for  even  moderate  values  of  vis¬ 
cosity  the  Voigt  model  will  act  as  a  rigid  element  in  its  im¬ 
mediate  response  to  high  magnitude,  short  changes  in  force. 

For  a  less  abrupt  rise  the  wave  velocity  and  the  rise  time 
approach  those  of  an  elastic  wave. 

Rise  and  Fall  Time.  The  spherical  pulses  in  the  Voigt  model 
show  a  marked  symmetry  in  their  first  positive  phase.  These 
are  similar  in  shape  to  the  pulses  in  GNOME  (salt)  and  AIR- 
VENT  (andesite),  (Figs.  3S  to  43)  and  are  also  symmetrical 
in  ti.eir  first  positive  phase.  The  pulses  for  HARD  HAT 
(granite)  and  LONG  SHOT  (andesite)  on  the  other  hand  have 
a  shorter  rise  time  and  a  longer  fall  time.  Exact  scaled 
comparisons  cannot  be  made  at  the  present  stage  of  develop¬ 
ments  with  the  Voigt  results  primarily  because  of  the  uncer¬ 
tainty  of  the  choice  of  a  radius  of  equivalent  cavity,  r, 
at  which  either  elastic  or  viscoelastic  behavior  could  be 
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assumed  to  become  effective.  The  Voigt  equation  is  linear, 
but  will  scale  only  if  the  viscosity  varies  as  the  scale  fac¬ 
tor.  This  condition  does  not  occur  for  small  and  large  events 
in  the  same  earth  material. 

It  would  appear  that  for  somewhat  porous  types  of  rock 
such  as  those  of  GNOME  and  AIRVENT  that  certain  properties 
of  the  .'ock  mass  cause  a  damping  out  of  the  higher  frequencies 
more  rapidly  than  in  granite  and  andesite.  Porosity,  water 
content,  and  geologic  structure,  as  well  as  elastic  proper¬ 
ties  of  intact  specimens  all  have  their  effect  in  varying 
but  unknown  degrees.  No  quantitative  coefficients  for  the 
damping  properties  of  large  rock  masses  for  high  order  tran¬ 
sients  have  been  found  in  the  literature,  however.  Hence, 
this  would  be  a  fruitful  area  for  further  research. 

The  value  of  c  for  the  wave  velocity  employed  in  the 
Voigt  wave  calculations  was  20,000  fps ,  or  a  value  equal 
to  that  of  a  dense  granite  with  large  shear  and  Young's 
moduli.  However,  the  shape  of  the  Voigt  pulses  corresponds 
more  nearly  to  that  of  low  moduli  materials  such  as  salt. 

Attenuation ■  Both  the  relative  magnitude  and  the  rate  of 
attenuation  of  the  normalized  particle  velocity  pulse  are 
affected  by  the  value  of  u0  (Fig.  34).  For  larger  values 
of  o o,  i.e.,  for  smaller  viscosity,  the  normalized  peak 
particle  velocity  becomes  larger  and  the  attenuation  is 
lower.  The  rate  of  attenuation  decreases  somewhat  with 
travel  distance,  which  is  not  the  case  for  actual  velocity 
pulses  in  salt,  granite  and  tuff.  However,  the  slope  of 
the  peak  velocity  curve  for  larger  values  of  w0  is  some¬ 
what  lower  than  that  of  observed  values  in  natural  rock, 
which  is  about  1.65  (Fig.  49). 

It  should  be  noted  that  the  calculated  peak  velocities 
are  for  a  wave  input  form  of  P(t)  ■  Poe"*5*,  which  has  an  in¬ 
finitely  steep  rise  in  pressure.  The  shape  of  the  input 
pulse  affects  the  attenuation  as  well  as  the  magnitude  of 
the  peak  velocity  as  shown  above.  At  the  radial  distance 
from  the  disturbance  at  which  either  the  elastic,  visco¬ 
elastic  or  similar  wave  equation  becomes  effective  (radius 
of  equivalent  cavity),  the  shape  of  the  pulse  crossing  this 
surface  is  also  affected  by  the  properties  of  the  rock  medium. 
That  is,  pressures  generated  by  explosive  sources  are  so  in¬ 
tense  that  there  may  be  vaporization  and  melting  (nuclear 
explosives)  as  well  as  pulverization,  fracturing  and  flow. 

The  limit  of  this  volume  is  described  by  the  distance  called 
radius  of  equivalent  cavity,  and  energy  and  pressure  losses 
will  be  determined  by  the  porosity,  strength  and  various 
properties  of  the  rock  which  contribute  to  its  ability 
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to  successfully  sustain  an  intense  pressure  pulse.  Thus, 
fractured  granite  apparently  can  sustain  a  pulse  with  sharper 
rise  time  than  porous  rocks  and  salt,  as  well  as  transmitting 
the  same  general  shape  of  pulse  beyond  the  fracture  zone. 
Further,  the  pulse  to  be  considered  for  possible  viscoelastic 
applications  would  have  a  more  gentle  rise  time  than  a  simple 
exponential.  Pulses  of  the  form  P(t)  =  P0  (e~at-e* ®t)  will 
be  investigated  in  a  continuation  of  the  study  of  viscoelastic 
waves . 

Conclusions .  The  Voigt  spherical  wave  generated  by  a  single 
decay  exponential  forcing  function  exhibits  some  characteris¬ 
tics  which  are  similar  to  those  of  pulses  observed  in  some 
kinds  of  natural  rock  masses.  The  Voigt  pulse  form  tends  to 
become  symmetrical  and  lengthens  somewhat  with  travel  distance. 
However,  the  velocities  of  the  frequency  components  of  the 
Voigt  pulse  are  frequency  dependent,  and  hence,  the  wave 
spreads  in  both  directions,  rather  than  lengthening  from  an 
arrival  time  described  by  the  elastic  velocity.  Observed 
pulses  in  salt  and  porous  andesite  likewise  tend  to  become 
symmetrical,  but  lengthen  more  than  the  Voigt  pulse  does  and 
only  in  one  direction.  The  attenuation  of  this  type  of  visco¬ 
elastic  wave  is  somewhat  greater  than  the  average  of  attenua¬ 
tion  rates  observed  in  rocks  in  which  contained  nuclear  explo¬ 
sions  have  been  carried  out,  but  is  comparable  for  tuff  and 
salt.  Waves  generated  by  pulses  with  longer  rise  times  will 
undoubtedly  more  nearly  approximate  waves  in  natural  rock. 

For  all  values  of  u0  investigated  and  for  the  single  ex¬ 
ponential  decay  pressure  function  the  arrival  times  of  Voigt 
waves  are  considerably  smaller  than  for  either  observed  or 
elastic  waves.  This  is  due  to  the  abrupt  rise  in  magnitude 
of  the  pressure  function  and  to  the  rigid  response  of  the 
viscoelastic  model  to  abrupt  changes  in  face. 

Thus,  while  a  change  in  the  pressure  function  will  make 
the  Voigt  wave  more  closely  approximate  natural  waves,  a  de¬ 
crease  in  the  viscosity  factor  causes  the  viscoelastic  pulse 
to  approach  the  character  of  an  elastic  pulse,  which  does 
not  lengthen  with  travel  distance.  It  appears,  therefore, 
that  while  for  certain  parameters  the  Voigt  mode’  may  more 
closely  approximate  observed  pulses  in  attenuation  and  pulse 
lengthening,  it  does  not  offer  sufficient  advantages,  except 
for  a  few  special  <.ascs,  to  justify  its  use  for  predicting 
wave  parameters  for  design  purposes  in  comparison  with  elastic 
waves  or  empirical  procedures. 

These  analyses  of  the  Voigt  model  indicate  that  a  con¬ 
figuration  other  than  one  incorporating  a  dashpot  clement  will 
be  required  to  repress',  t  rock  materials. 
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TABLE  I  *  Rock  Properties 


Physical  Property  Data 


LONG  SHOT  -  andesite  (Ref. 

12) 

Compressive  strength: 

2260  - 

5260  psi 

Specific  gravity: 

2.27  - 

2.36 

Ultrasonic  pulse 

velocity : 

6100  - 

10,645  fps 

HARD  HAT  -  granite  (Ref.  13) 


Compressive  strength:  10,835  psi 
Tensile  lit  strength:  1,915  psi 
Bulk  de  ..xty:  2.69 

Ultrasonic  velocity 

(20  kc) :  19,450  fps 

Modulus  of  elasticity ,E:  11.3  x  lO^psi 

Poisson's  ratio:  0.20 


GNOME  -  salt  (Ref.  14) 


Compress ional  wave 
velocity : 

Shear  wave  velocity: 
Density  (dry  ,  bulk) : 
Poisson's  ratio: 
Young’s  modulus: 
Porosity  (1  volume): 
Water  content  (1  wt) : 


4.08  km/sec 
2.88  km/sec 
2.13  -  2.46 
0.28  -  0.3’ 

2.41  x  101 *  dynes/cm^ 
0.76  -  5.H 
11 


AIRVENT  -  data  not  available 
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Cttffiiaa  FLOW  .SHEET 


Fig.  1.  Generalized  Computer  Flow  Sheet 
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Fig.  5.  Normalized  strain  e'(x,t)  ■  e(x,t) 
small  x  for6(t)  forcing  function. 
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spherical  Voigt  wave,  r0  ■  50  ft.,  (d<>  “  600  and  c  “  20,000  ft/sec 
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I _ I _ I _ I _ l 

0  10  20 

time-  ms 

fl9-  33*  Normalized  particle  velocity  v'(r,t)  *  v(r,t)  :  ■■**■'■* 

3uw.»  n 

for  P(t)  -  P.e  for  spherical  Voigt  wave,  r.  -  50  ft., 
u>.  -  3000,  S  -  150,  and  c  ■  20,000  ft/sec.  Notation  t 
indicates  arrival  time  of  elastic  wave. 
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*nd  r  »  150. 


-  ft/sec 


0 


0.2 


0.4 


0.6 


time  -  sec 


Fig.  36.  Particle  velocity  pulse,  Project  A1RVENT,  gage  60u. 
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37.  Particle  velocity  pulse,  Project  A1RVENT,  gage  70u. 
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time  •  sec 

39-  Particle  velocity  pulse,  Project  AIRVENT,  gage  130-u. 
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m/s 


Fig.  4).  Particle  velocity  pulse,  Project  GNOME,  range  229  meters,  6-AH. 
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0  0.2  04 

time  -sec 


Fig.  43.  Particle  velocity  pulse  in  salt.  Project  GNOME , 
range  477  meters,  21-A1H. 
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ft/sec 


time  -  sec 

Particle  velocity  pulse  in  andesite,  Project  iONGSHOT, 
80  KT,  h  irizontal  range  600  ft.,  gage  depth  100  ft., 
s’ant  range  2281  ft.,  gage  ''H,  EH-?. 


1.0 


10 


100 


range,  100  ft/kt 


Fig.  49.  Composite  correlation  of  magnitude  (peak)  of  radial 
velocity  pulse  in  salt,  granite  and  tuff,  scaled  to 
1  KT  (Ref.  15). 
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APPENDIX  A 


WAVE  PARAMETER  FORMULAS 


Summation  Symbols 
Plane  wave  -  single  summation 


(t’P)-  V  X^,(n*z)/2  (n+?.-l)/2 

z  L  n!  -n-z 

n=o 


Plane  wave  -  double  summation 


II 


(s,P) 

y 


CO 


I 


03 


I 


expC-Xv's  j 
sm+n/2+y 


I  lf'P)=  I  I  bm  2m^n+z)/2  Tm'f (n+z*l)/^D  /  X  \  Exp 

L  Lz  L  L  m  n!  -n-2m-z  V/tt/ 

m=o  n=o  V*'/ 


Spherical  wave 

(S»s)^  "  “  gXB.C-Rjs) 

L  Lv  L  L  m  n!  (m+nl/2+y 

1  m=o  n=o  sv  "  ’ 


I  I 


(t.S) 

z 


03  00 


I  I 


b 


m 


Rn  7(ni+n+z)/2 
n! 


T(m+n+z-l)/2D 

-n-m-z 


(s,P)  =  s  variable,  plane  wave  (s,S)  =  s  variable,  spherical  wave 

(t,P)  =  t  variable,  plane  wave  (t.S)  =  t  variable,  spherical  wave 

On  the  left  side  of  the  above  equations  the  first  letter  in 
the  superscript  represents  the  s-plane  or  the  t-plane,  and 
the  second  a  plane  or  spherical  wave.  The  subscripts  y  and 
z  a  e,  respectively,  the  numerical  values  of  the  transform 
variable  exponent  other  than  the  summation  indices,  and  the 
subscript  of  the  Weber  function  other  than  the  summation  in 
dices.  In  all  cases  z  =  2y-l,  EXP  *  exp(-T-X2/8T)  for  the 
plane  wave,  and  EXP  «*  exp( -T-R2/8T)  for  the  spherical  wave. 


83 


^  is  the  transform  solution,  (2)  the  transform 
solution?”  nornalized'  (3)  the  expanded  transform,  and  (4)  the  inverse 

Plane  Wave:  P(t)  =  P'6(t) 

Particle  Velocity 


(1)  v(x,t)  *- 


*-  P'  t.»xp[-xs/c(s/Mo>i)^] 


pc 


(s/oio  +  l)  ' 


(2)  (1/ui,)  vCx.t/m.Je*  '=  P-exp(-X/r,X^) 

pc  7s 


(3) *  (l/uo)  V(x,t/aj„)et  *=  —  j[S,P') 

pc  Li 

(4) *  v(x, t)  =  ^2.  y(t,P) 

pc/ir  Lo 

Displacement 

Cl)  u(x,t)  P'  exp[-xs/c(s/u0+l)*] 
pc  s(s/u),+l)i 

(2)  (l/ooe)  uCx.t/woje*  *“  —  exp(-x7r»X/7s) 

pcwo (s-1)  7i" 

(3) *  ufx.t/Me*  =  tL  l  I  S,P) 

pc  ^  2.3/2 

p.  (t.P) 

(4) *  u(x, t )  =  P  l  I 

pcni  L  L2 


*  In  summations  all  b  =  l 
m 
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Strain 


(1)  e(x.t)  H 


P'  exp[-xs/c(s/u)0+l)^] 
pc2  (s/<Do  +  l) 


(2)  (1/uo)  eCx.t/woje1  =. 


P1  exp(-X/s+X//s) 
pc2  s 


(3)*  (1/oj.)  etx.t/ajoJe1  =. 


p.  (s,P) 


(4)*  e(x,t) 


P'too 

pc2/jT 


(t.P) 

i 


Stress 

(1)  a (x, t)  H  -  P'exp[-xs/c(s/a)0+l)^] 

(2)  (1/ojo)  (x.t/aioje1'  H  -  P'expf-X^X/Zs') 

r  (s,P) 

(3) *  (1/wo) (x,t/wo)er  =  -  P*  lQ 

r(t.H 

(4) *  a  (x,  t)  -  -  P’ooo  l_ 


**  In  summations  all  b  =  1 

m 
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Pljine  Wave:  P(t)  =  P"l(t) 

Particle  Velocity 

(1)  v(x,t)  P”  expt-xs/cCsWl^] 
PC  (s/w.+l)* 


(2)  (1/ Wo)  vtx.t/a^e1  =  El.exp(-X/s*X/.£). 

PCWo  (s  - 1)  /s' 


(3)*  Cy  £<s-l>) 

PC  “  LZ/2 


nil  (t#P) 

(4)*  v(x,t)  =  —~l  l 
pc/n  2 


Displacement 

(1)  u(x,t)  Up"  exp[-xs/c(s/a,o+l)*l 

pc  S2/s/uio  +  l 

(2)  u (x# t/ u)o) el  v—.  expC.x/5Vx/^) 

pCaio  (S-l)2^ 

f  p,f  (s  t  P) 

(3) **  uCx.t/^.Je1  - -  l  I 

PCWo  *■  *•  3/2 


(4) **  u(x,t)  « 


P^ _ 

PCw«/iT 


*  In  summations  all  b  «  1 

n 


Summation  on  m  is  from  1 
and  in  summations  b  *  m 


for  this  case  only 


Strain 


(1)  e (x, t)  =  -  p-  exp[-xs/c(sA>.+l)*1 
pc2  S(s/ui„  +  l) 


(2)  (1/w.)  etx.t/u.Je*  *-  -  P^ex^-X^X/ZT) 

pc2u0  S(s-l) 

(3) *  (1/u.)  eCx.t/w.Je1  £1L  V  y(S'P) 

P C*  L  L2 

pit  (t ,  P) 

(4) *  e(x,t)  =  -  I  l 

PC2 /it  l  l  3 


Stress 


(1)  o fx.tl  •-  -  p"  exp[-xs/c(s/w.»l)*l 

s 

(2)  (1/w.)  olx.t/w.Jc1  •=-.  -  -P"  -  exp(-x/5~+X//T) 

w«(s-l) 


(3)*  JCx.t/iuole1  H 


P"  l  I 


(s.P) 

l 


„  (t.P) 

(4)*  0(x,t)  -  -  P"  l  l 


*  In  summations  all  b  *  l 

m 
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Plane  Wave:  P(t)  =  P0e 
Particle  Velocity 


■et 


(1)  ,(x,t)  St 

PC  0\  <c  >.11* 


(s+6)  (s/'j)o+l)‘ 

.  t  .  Po  expl-X^s+X/Zs) 

(2)  (1/u.o)  v(x,t/u»e)e  pcWo  (s.a)Vs'_ 

where  a  *  (1-6/uo) 

t  P  (S*P) 

(3) *  v(x,t/a)o)e  ^  l  I3/2 


(«)•  v(,.t> 


(t.P) 


Displacement 

(1)  »«,.)  p. 

pc  s(s*6)  (s/u».*l) 

t  P.  exp(-X»^*X //s’) 


(2)  u(x,t  /usoje 


o cuio  (s-1)  (s-ci)  »s 
(S.P) 


(3)**u(x,t/ii)«  )c  * 

P.  T  T 


t .  p.  V  y  ■ 

CJ.  i 

lt,P) 


(4)**u(x,t) 


P  Cw« 11  n 


•  In  summation  b 

a 


a 

v.l 


•  • 


In  summation  b  obtained  fro»  expansion  of  l/(s-l)(s-a) 


Strain 


(i)  s(x,t)  h 


P0  e  xp  [-x3/c(3/aio+l)^} 
pc2  (s+8)  (s/aio+1) 


(2)  e  (x,  t/a»o)et 


P8  exp  (~x/s+X//s~) 
pc2  s(s-a) 


(3)*  e(x,t/wo)et 


l  l 


(s,P) 


(A)*  e (x, t) 


(t,P) 

3 


Stress 


(1)  a(x,t) 


P„  exp  [-xs/c(s/u)o+l)^] 
(s+ci) 


(2)  c(x,t/..)ec \  - 

(s-a) 

t  (s,P) 

(  3)  *  P(x,t/w.)er  \  -  P,  )'  l 


(A)*  J ( x , t  > 


(t.P) 


*  In  suncuUicns  b 

o 


m 

i 
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Spherical  Wave:  P(t)  *  P'6(t) 


Displacement 


(1)  ucr.t)  [jl  ♦ _ ? _ 1 

u(s/ti)«*l)  [B]  rc(s/wo*l)H 


(i)*  u(r.t/u.)etH  ^£- 

3  MU* 


,  (S.S)  (s,S)  (s,S) 

V  n2  *u-n3/2-^ns/2 

2  rc  3/2  rc  5/2 


p  •  * 

(4)*  u(r.t)  - 

3uw«  r* 


,  U.S)  (t,S)  (t.S) 

*it.yy  i±s.  y  y 

^ ’  rc  L  L2  rc  L  Lu 


-3 


*  In  summations  all  obtained  from  expansion  of  1 / ( A ] 
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Particle  Velocit 


(1)  v(x,t)  H  p'r«s  exp[-s(r-r.)/c(s/h>,»l)*]  fl  ^  s 

u(s/u«+l)  [B]  Jr2  rc(s/o)0+l)i 

(2)  p,r»g2(s-D  exp(-R^7+R/,£)  fj_  ^  u.(s-l)  1 

3g[A]  '  2  ~~tc7s~  J 


In  summations  all  obtained  fro*  expansion  of  1/{A) 
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Strain 


I-- 


m  ,  >  •_  P'r0  exp[-s(r-ro)/c(s/a)o+ii5} 

p(s/li)e  +  l)  [BJ 


O  ‘ 


f*— -  ♦  - — - r  ♦ - *- - ] 

[_r3  r2c(s/tu0+l)  ‘  rc2(s/u«+l)J 


(2)  Efr.t/w.je1  *=  -  L^oclexLL:^l/2^. 

3yaio  [A] 


[2_+  20)0  i 

|_r3  r2ci 


S-l)  +  b)oZ(S-l] 


t.  P'rc2  (2  (s*S>  2„  Fp  p(s*S)  p  Js>S)1 

(3) *  e(r,t/u0)e  ^  l  l2  *  jl  ~  l  I$/2  j 

„  2  r  (S,s)  p  (s,S)  (s,S)l) 

+  -2^2  *ns  JJ 

P'r  C2  f 2  r  r(t'S)  2m  fp  p(t>S)  p  r(t.S)"l 

(4) ‘  eM  =  -^p”3  *^[n2  -n4  j 


2  T  it.S)  (t,s)  (t.s)!) 

2*  II,  -2^3  .II,  Jj 


In  summations  all  obtained  from  expansion  of  1/ [A] 


Stres? 


(1)  o(r,t)  n  -  — 


?*ro  exp[-s(r-r«)/c(s/fc>o->-l)*] 
5  [31 


pL  ♦ _ 1! _ T  ♦  — il2- _ 1 

|_r-  r2c(s/a)o+i}a  re2 (s/u0+\)  J 


(2)  a  (r,t)et  H  -  g’^s.e^-Rvgig^gS. 
r  3o;o  [A] 


f-i  .  4w0 (s-1)  .  3<uo2 (s- 1)21 

F  XX 


(!)■  *  (,.«).« •-  -  £n£  ti.  I  f’\  if  ft  y(s,S)- 1  y(s,s  1 

r  3^,,  U3  L  Ll  r2c  [  “V2  3/2  J 

*  g#  [i  ci,S)-  2 1  i!s'S)*  i  i!s'S)]} 

T.,  2  r„  „(t>S)  _  _l'£»S) 

*  S3- 1 1,  - 2  r  i,  *  a,  \f 


I 


'  In  summations  all  bm  calculated  from  expansion  of  1/ [A] 


1  ><rK,^,nv.i'^V:S'»i>ivij'.;ri,Vv-l,Vv  •• . 


Spherical  Wave:  P(t)  =  P"  l(t) 


Displacement 


(1)  u(r,t)  -  P"r<>  exp[-s(r-ro)/c(s/ue+l)*] 

psCs/'uo+l)  [B] 


1  s 


r2  rc(s/w0+l): 


(2) 


ufr.t/aioje1  H  -"ZsS*.  exp(;-R^R/^)  [l_  + 

A]  r2 


3po)02(s-l)  [A]  j  r2  rc/s" 


(3)-  P"r"c2  <•  '(S,S>-  “• 


.  (s  ,S)  (s,S) 

3pu)o2  I  r2  “  “3  +  „  P  5/2  '  ^  7/ 2 


(4)*  u(r.t)  =  ~1’£—  l  ft,S)+  £2- 

3wfc:«.  i/ir  I  r2  5  rc 


l\  Ct.S)  (t,S)’ 

1 1  - 1  L 

4  6 


*  In  summations  bm  calculated  from  expansion  of  l/(s-l)[A] 
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Particle  Velocity 


(1)  v(r,t)  « 


♦.  p"to  exp[-s(r-r<,)/c(s/aie*l)^]  fl 


y(s/u0+l)  [BJ 


L?‘ 


rc(s/u0+l) 


T 


(2)  v(r>t/a..)et  exp(-R^V/^)  fl  +  Ms-l)l 

3mu!o  [A]  It2  tc/s 


<•  pti-  .2  (i  (S»S)  (S, S)  (S 

(3)*  v(r,t/a)o)et  =,  )L-  V  £  +  <Sl  H  .TV 

3ym.  )r2  L  L2  rc  L  Ll/2  L  lS/2 


;"]! 


P.._  .2  ( ,  „  „(t,s)  r  ct,S)  (t.srn 

w’ 5=77(7”.  'Z[n,  -Ik  ]( 


*  In  summations  bm  calculated  from  expansion  of  1/[A] 
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* 


Strain 


(1)  c  (r  1 1  -  P"r°  exp[-s(r-ro)/c(s/a>.-»l)*] 

uCs/aio+1)  [B] 

\i - 22 —  — ii — 1 

Lr3  r2c(s/u0+l)S  rc2(s/a)o+l)J 

(2)  efr.t/o.oje1  =  - 

3mioo2  (s-1)  [A] 


2  +  2u)0  (s-1)  Wo 

.  r2  r2c/s" 


(3)*  £(r,t/ui0)et  — ~ 

3uu)o  i 

2  r„  Js>s) 


2(S-1)21 

rc2s  J 

(?  r  T.^S»S^  JS>S)  „  n(S,S)*| 


Wo 

rc2 


(4)*  cCr.t)  = 


Mo* 

j *r>2 


p  p,  -  -  p  (s,S)  (s.S)-li 

’i,  '2”3  ]} 


P"roC2  (2  y  y(t'S)  2wo  fy  y(t>s)  y  y(t’S)l 

7^777  p  1  S  +  7^  [lK  ""  6  J 


5 

r  (t,S)  (t,s)  ,  (t.S)  \ 

[”,  -2”s  J| 


*  In  summations  all  calculated  from  expansion  of  1 / [A ] 


96 


Stress 


(1)  ci(r  t)  =  -  P"  exp[-s(r-r8)/c(s/oie-t-l)jr] 

■’  s  [B] 

f4  +  4s  +  3s2  ] 

Lr7  r2c(s/o)0+l)i  +  rc^(s/oic+l)  J 

(2)  o(r,t/u0)et  *=.  -  ^>-C--S  ^xp(-R^R//7)  x 

3w02  (s-1)  [A] 


i-  +  4Ms-P  + 

r3  r2c/s" 


3gj02  (s-I)  21 
rc2s  -1 


*  In  summations  all  calculated  from  expansion  of  1/ (A] 
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Spherical  Wave:  P(t)  »  P«e 


-St 


Displacement 

(1)  u(r,t)  P°r«  •«Pt-»(i^g>)/c(B/m.-l-l)1*]  T~i_  + _ g_ 

w(s+8)  (s/aio+l)  [B] 


i_  +  S 
r2  rcCs/uo+l)1* 


(2)  u(r,t/ui.)et  H  b**?  exp(-R/^-R//j) 

3pw02  (s-a)  [A] 


1  .  UJo(3-l) 

r2  rc/s 


a  -  (1-8/ui.) 


(3)*  utr.t/auje*  H  j ±r  Y  [(s*S)+  SO. 

3uaio‘  (r2  L  S  rc 


l  [(S,S>-  l  l 3 

*  S/2  *  S/ 


:"]) 


(4)*  u(r,t)  - 


3uu>©2«'/iT  lr2 


ii  .  Sfc>s) 

f(?n  ♦ 


rj  i,<t,s)- i  i<t,s> 


*  In  summations  calculated  from  expansion  of  1/ (s-a) (A ] 
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.  Particle  Veloclt 


(1)  v(r,t) 


L 

P pros  exp[-3(r-r«)/c(s/a)a+l)  ] 
u(s+6)  (s/aio+l)  [B] 


1_  +  8 
r2  rc(s/w«+l) 


(2)  v(r,t/uio)et  h 

3uu0  (s-a)  [A] 


i_  +  ui.(s-l) 
r2  rc»^ 


►  d  JL  (8,S)  (s  ,S) 

»>*  &  n2  -n, 


rc 


(s,S)  (s,S)  <s,S) 

1 2,/,  - 2?  *  J  I,„ 


“>*  v<r-t)  •  wE  {ir 


(t,S'  (t,S) 

1 1, 


+  2* 

rc 


(t,S)  (t,S)  (t,S) 

n,  -  *n.  ♦  i l 
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Strain 


(1)  c(r  t)  P«re  exp[-s(r-ro)/c(s/ui0»l)*] 

u(s+S)  (s/uo^l)  [B] 


3  ~ 

r3  r 


2s  s' 

4- 


2c(s/(Do-*’l)^  rc2(s/a)o+l) 


(2)  cfr.t/Me1  *-  -  fl_  +  2o^(s_-lj_  +  y.2_(5-l)21 

3y(s-a)  [A]  |_r3  r2c</s"  rc2s  J 


»t  ‘m.  Ppr»c2 


{ 2  r  r(s.S)  .  r  (s,S)  (s,S)-| 

(3)*  e(r,t/u,.)e'- ’■?-  <MI  ♦#L  ll  -II 

3yu.2  |r3  L  S  r2c  S/2  **  S/2 


2  r  (S>s)  (*,s)  (s,s) 

Mu,  -ui,  -u. 


]) 


(4)*  e (r,t)  = 


o  2  {->  (t.s)  .  r  (t,S)  (t,S)-| 

'f^r-  Ml  l  l  'll 

3inii.Vi  {r3  5  r^c  L  4  6 


c»a 


2Tr  r^t,S^  r  r(t*S)  rr(t*S)l) 

[u,  -«n,  *n;  J) 


rc 


In  summations  b  calculated  from  expansion  of  l/(s-a)  [A] 

UR 


Stress 


(1)  o(r  t)  -  -  p«r«  exp[-s(r-ra)/c(s/ui.t-l) 

(s+8)  [B] 


[i -♦ _ «* _ _  ♦ _ 151 _ 1 

j_r3  r2c(s/«Jo+l)i  rc2 (s/wo+l) J 


(2)  c(r, expl-B^R/^ 

3«o2  (s-a)  [A] 


+  Au)«  (s  - 1)  +  3u)q2(s-1)21 
_r3  t^o/s  rc2s  J 


t.  ,t._  P.r.C2  |  4  r  r(S'S)  4W.  fy  r(S,S)  t  f$' 

(3)  a(r.t/«.)e  •  ~  fa  l  l2  ^  [^3/2  '  H5/2  J 


*  In  suiBMtions  calculated  fro«  expansion  of  l/(s-a)[A] 
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Plane  .  mi  spherical  waves  in  a  Voigt  medium  were  mve iga ‘ni  to  .empar* 
the  calculated  wave  tor  ms  wi'h  observe*!  waves  generated  bv  arg»  contained  Hi.  and 
NE  explosions.  Interest  is  center*  d  on  wave  terms  m  \vh.u  is  usually  v  omnirrn!  to 
be  the  I'iastis  region  around  an  explosion.  Plane  waves  do  not  apply  at  this  distance 
hr  c. iu«te  of  geometry.  The  plane  Voigt  wave  equr.tmn  has  beer,  previously  solved  for 
particle  vc’  v  1 1  v ,  stress  and  strain  tor  a  unit  impulse  forcing  furu  lion.  However, 
solution*  lor  the  dif  pUc  r tr-nf  for  a  unit  impulse  and  for  the  four  wave  parameters 
«t»r  a  unit  step  and  a  dr*  av  exponential  involve  multipliers  in  t.ie  operational  form  fur 
which  no  transform  pairs  have  been  published.  A  method  of  solution  i*  pr»s  ntrd 
uliii  h  utilises  a  Heaviside  expansion  of  the  multiplier*  in  the  transform  plane  which 
result*  *n  products  of  two  mimitr  series  whu  h  mu\  be  inverted  tern:  bv  tern?.  These 
may  hr  Ivirlher  resolved  a#  Single  series  with  polynomial  inefficient*  C  r  purposes 
t*s  i  o.nputation,  •*.  similar  method  of  solution  of  *«ugt  spherical  wave*  w.»*  found  for 
ani!  itiip  il*c,  unit  step  ar.d  dr  a\  exponential  Ion  ir^  functions  f«*r  displ  a  emrnt ,  par- 
tide  velocity*  and  radial  stress.  Appropriate  recursion  tormula*  make  them 

readily  adaptable  to  computer  evaluation.  O**  illation*  **u'  ir  for  a  sphere  al  Wave 
whereas  for  a  plane  wave  they  do  not.  (  a  U  >1  ation .«  were  performed  for  parti*  ie  ve¬ 
locity  ;<*r  three  value*  of  the  Voigt  viscoelastic  parameter  „  and  »  c  -nt'a r « son*  made 
usth  pulse  f.ifm*  for  waves  in  granite,  tuff  an-  salt. 
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Please  make  the  following  revisions: 

1.  Insert  on  page  33: 

"AIRVENT*  -  playa 

Density  -  wet  1.44-1.85 

Density  -  dry  1.15-1.67 

Porosity  -  estimated  5070 

Hugoniot  Data 

Pq  =  1.41-1.47  gm/cnP 

V  =  2.58-5.24  mm/iisec 

u  =  1.04-3.54  mm^sec 

P  =  30-271  kb 

v/v  =  0.302-640" 

'  o 

2.  Page  28,  line  10  -  Delete  "AIRVENT" 

Insert  after  line  10  -  "AIRVENT  was  a  40,000  lb,  TNT  event 

detonated  at  17  ft.  depth  in  desert  playa." 

3.  Page  29,  line  26  -  Delete  "No  data  were  available  for  AIRVENT." 

4.  Page  30,  line  37  -  Substitute  "playa"  for  "andesite." 

*  Kintzinger,  P.  R.,  AIRVENT  Phase  I,  Project  I,  Earth  Particle  Motion. 
SC-RR-64-549*  Sandia  Corporation,  October  1964. 


